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^ ■ ABSTRACT 

\ We present the optical identification of mid-IR and radio sources detected in the 

-«.^ , European Large Area ISO Survey {ELAIS) areas Nl and N2. Using the r' band optical 

■ data from the Wide Field Survey we apply a likelihood ratio method to search for the 

counterparts of the 1056 and 691 sources detected at 15//m and 1.4 GHz respectively, 

Q . down to flux limits of 5*15 — 0.5 mJy and 51.4 ghz — 0.135 mJy. We find that '^92% of 

$H ' the 15/Ltm ELAIS sources have an optical counterpart down to the magnitude limit of 

\ the optical data, r'=24. All mid-IR sources with fluxes ^is ^ 3 mJy have an optical 

■ counterpart. The magnitude distribution of the sources shows a well defined peak 

at relatively bright magnitudes r''^18. About 20% of the identified sources show a 
point-like morphology; its magnitude distribution has a peak at fainter magnitudes 

' than those of galaxies. The mid-IR-to-optical and radio-to-optical flux diagrams are 

Jh \ presented and discussed in terms of actual galaxy models. Objects with mid-IR-to- 

^ optical fluxes larger than 1000 are found that can only be explained as highly obscured 

star forming galaxies or AGNs. Blank flelds being 8% of the 15/im sample have even 
larger ratios suggesting that they may be associated with higher redshift and higher 
obscured objects. 

Key words: galaxies: infrared: galaxies - galaxies: evolution - star: formation - 
galaxies: starburst - cosmology: observations 



1 INTRODUCTION 

The Infrared Space Observatory (ISO, iKessler et al]|l996l) 
was the second infrared space mission, providing a great 
improvement in sensitivity over the IRA S mission. The Eu - 
ropean Large- Area ISO survey {ELAIS, lOliver et aLlboOOH 
was the largest Open Time programme on ISO. This project 
surveyed 12 square degrees, divided into four main fields, 
three in the north (Nl, N2, N3) and one in the south (SI). 
The main survey bands used were 6.7, 15, 9 and 170 /im. 
The ISO CAM camera dCesarskv et al.l 119961) was used for 
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the s horter wavelengths while the ISOPHOT jLemke et alJ 
Il996h camera was used for the longer ones. 

Optical imaging is essential to study the properties of 
the sources detected. Due to the large errors ellipses of the 
mid-IR detections, typically several seconds of arc, it is nec- 
essary to carry out a detailed process of identification. Since 
more than one optical source can be inside those ellipses, 
a method which provides the likelihood of each counterpart 
to be the true association is needed. The identification does 
not only provides us the optical properties of the mid-IR 
sources but also allows us to further proceed with foUowup 
observations of interesting sources. 

This paper presents the optical identification of the mid- 
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Figure 1. Schematic layout of the Wide Field Camera. All the 
measures are approximated. 



IR and radio sources in the Nl and N2 areas. Section |21 
presents a summary of the optical observations carried out 
in these areas as well as the reduction steps and products. 
Section |3| describes the mid-IR and radio catalogues used. 
Section 0| discusses the actual procedure to determine the 
optical counterparts of the sources, while sections |K| and |S| 
describe the optical properties of the sources. 



2 THE OPTICAL CATALOGUES 

In order to identify the mid-IR sources with optical ob- 
iects we use the data from the Wide Field Survey (WFS, 
iMcMahon et al.ll200ll) . This survey has been carried out 
using the Wide Field Camera (WFC) on the 2.5 m Isaac 
Newton Telescope (INT) on the Observatorio del Roque 
de Los Muchachos (La Palma). The WFC is formed by 4 
4kx2k CCDs. The arrays have 13.5 /im pixels correspond- 
ing to 0.33" /pixel at the telescope prime focus and each one 
covers an area on sky of 22.8x 11.4 arcmin. The total sky 
coverage per exposure for the array is therefore 0.29 square 
degrees. 

Figure shows an schematic layout of the CCD detec- 
tors in the WFC. Gaps between detectors are typically 20". 
Chip 3 is slightly vignetted in one corner. Optical observa- 
tions described below are carried out considering the WFC 
as a 3-CCD camera, allowing for a 10% overlap between 
adjacent pointings for photometric purposes. Therefore the 
area not observed in each pointing due to the chip gaps is 
about 12 square arcmin. However, the spatial source density 
of ELAIS and radio sources is not high enough to make this 
effect important for the process of identification. 

The WFS surveyed ~200 square degrees in different well 
known regions of sky with data at other wavelengths. The 
ELAIS regions Nl and N2 were also included. The survey 
consists of single 600 s exposures in five bands: U, g', r', i' 
and Z (see figure|2l to magnitude limits of: 23.4, 24.9, 24.0, 
23.2, 21.9 respectively (Vega, 5(t for a point-like object), 
i.e., about 1 magnitude deeper than the Sloan Digital Sky 
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Figure 2. Transmission curves of the U, g', r', i' and Z WFC 
filters (dotted line). Dashed line shows the quantum efficiency of 
the detector, where solid line shows to total throughput of the 
system. 



Survey (SDSS; lYork et al.ll2000h . A total of 108 pointings 
were done in Nl and N2, covering a total area of 18 square 
degrees. Typical seeing is about 1.0-1.2". The data are pro- 
cessed by the Cambridge Astronomical Survey Unit (CASU) 
as described in Irwin & Lewis (2000) and we provide here a 
short description of the reduction steps. The data are first 
debiassed (full 2D bias removal is necessary) . Bad pixels and 
columns are then fiagged and recorded in confidence maps, 
which are used during catalogue generation. The CCDs are 
found to have significant non linearities so a correction using 
look-up-tables is then applied to all data. Flatfield images 
in each band are constructed by combining several sky flats 
obtained in bright sky conditions during the twilight. Expo- 
sures obtained in the i' and Z bands show a significant level 
of fringing (±2% and ±6% of sky respectively). In order to 
remove this effect, master fringe frames are created by com- 
bining all the science exposures for each band. These fringe 
frames are then subtracted from the object exposures. Af- 
ter this removal, the fringing level is reduced to ±0.2% and 
±0.4% of sky in the i' and Z bands respectively. Finally an 
astrometric solution starts with a rough WCS based on the 
known telescope and camera geometry and is the progres- 
sively refined using the Guide Star Catalogue for a first pass 
and the APM or PMM catalogues for a final pass. The WFC 
field distortion is modelled using a zenithal equidistant pro- 
jection (ZPN; iGreisen fc Calabrettall2002^ . The resuhing 
internal astrometric precision is better than 100 mas over 
the whole WFC array (based on intercomparison of overlap 
regions). Global systematics are limited by the precision of 
the APM and PMM astrometric catalogue systems and are 
at the level of 300 mas. The object detection is performed 
in each band separately using a standard APM-style object 
detection and parametrisation algorithm. Standard aperture 
fiuxes are measured in a set of apertures of radius r/2, r, 
\/2»", 2r, 2\f2r where r = 3.5 pixels and an automatic 
aperture correction (based on the average curve-of-growth 
for stellar images) is applied to all detected objects. 

Photometric calibration is done using series of Lan- 
doldt standard stars (Landoldt 1992) with photometry in 
the SDSS system. For each night a zero point in each fil- 
ter is derived. For photometric nights the calibration over 
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the whole mosaic has an accuracy of 1-2%. During non- 
photometric nights, in otherwise acceptable observing con- 
ditions, we find that the derived zeropoint systematic errors 
can be up to 10% or more. Although the pipeline usually suc- 
cessfully flags such nights as non-photometric it still leaves 
open the problem of what to do about tracking the varying 
extinction during these nights. 

All calibration is by default corrected for the mean at- 
mospheric extinction at La Palma during pipeline process- 
ing (0.46 in U, 0.19 in g', 0.09 in r' and 0.05 in i' and Z). 
Since adjacent camera pointings overlap by several square 
arcminutes sources in these overlapping regions can be used 
as magnitude comparison points. However, to use overlap- 
ping and hence in general, different CCDs, requires that any 
repeatable systematics due to, for example, slight differences 
in the colour equations for each CCD, are first corrected for. 
Correcting for these is a three stage process. First the twi- 
light flatfields are used to gain-correct each CCD onto a 
common system. However, since the twilight sky is signif- 
icantly bluer than most astronomical objects, a secondary 
correction is made using the measured dark sky levels in each 
CCD for each filter to provide a correction more appropri- 
ate for the majority astronomical object. These corrections, 
unsurprisingly, are negligible for passbands on the flat part 
of the generic CCD response curves such as g' and r', and 
amount to 1-2% for the i' and z' passbands. The measured 
dark sky values for the U-band were also consistent with 
zero correction though with less accuracy due to the low 
sky levels in the U-band images. Finally any residual offsets 
between the CCDs are checked for each survey filter using 
the mean offset between adjacent pointings on photometric 
survey nights. The only filters requiring significant adjust- 
ments to the individual CCD zero points at this stage are 
CCDS for U (-3%) and CCDl for z' (-^3%). 

Data from non-photometric nights can now be cali- 
brated in one of two ways: the overlap regions between 
pointings can be used to directly tie in all the frames 
onto a common system, with extra weighting given to data 
taken from photometric nights; the stellar locus in vari- 
ous 2-colour diagrams (in regions of low unchanging ex- 
tinction like these) can be used to compare colours from 
all the passbands by cross-correlating the loci between any 
of the pointings. The colour-colour loci cross-correlations 
(actually made from a smoothed Hess-like version of the 
diagrams) give results accurate to better than -|-/-3% and 
can be used in conjunction with the overlap results to give 
overall photometry for the survey to the level of «2%. 
The final products include astrometrically calibrated im- 
ages as well as morphologically classified merged multicolour 
catalogues, publically available from the WFS web page 
(http : //www. ast . cam. ac .uk/^wfcsur /index . php). 

For the purpose of this work, we have also carried out 
the object detection in the r' band images using SExtrac- 
tor iBertin fc Arnorit3ll99£ri . As well as providing a inde- 
pendent check for the extraction method used above, it gives 
us an extense set of parameters. Aperture magnitudes from 
SExtractor are found to be in agreement with the ones ob- 
tained from the WFS pipeline. Therefore, together with the 
WFS aperture magnitudes we use MAG_BEST as a measure- 
ment of the total magnitude in the r'-band. The galaxy- 
star classification given by CLASS_STAR is also used to select 
point-like objects (defined to be CLASS.STAR > 0.8). 
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Figure 3. Comparison between WFS and SDSS r' band magni- 
tudes for point-like objects. A correction Am = 0.09 has been 
applied to account for the sUght filter and system diflferences. 



As an additional test of our photometric calibra- 
tion we have correlated the WFS catalogues with those 
from the Sloan Digit al Sky Survey First Data Release 
iAbazaiian et all 12003^ . Figure |3] shows the result of this 
correlation for the r' band magnitude. Once accounted for 
the small correction between both filters (Am = 0.09) the 
agreement is within 0.04 magnitudes (note also that a sec- 
ond order correction due to different object spectral energy 
distributions is not removed, so the accuracy is better than 
the quoted 0.04). 



3 THE MID-IR AND RADIO CATALOGUES 

The final analys is of the 15/xm data using the Lari method 
l|Lari et alJl2001^ for the ELAIS northern fields has recently 
been completed (Maccari et al. 2004, in prep). The two 
main northern fields Nl and N2 are centred at 16''10™0T' 
4-54°30'36" and 16''36'"58^ -F41°15'43" respectively. We 
have obtained a sample of 1056 sources (490 in Nl and 566 
in N2) to a 5(t fiux limit of 0.45 mjy. This catalogue includes 
sources detected in deeper observations (in particular, the 
central 40' x40' N2 area has been observed three times) so 
this flux limit is not homogenous over the whole survey area. 

As part of the multiwa velenth follow-up o bservations 
carried out in these regions ICiliegi et alJ 1^2^ have con- 
ducted a survey at 20 cm using the VLA in its C conflgura- 
tion, covering 4.22 sq. deg. in the Nl, N2 and N3 areas. They 
detect a total of 867 sources (362 in Nl, 329 in N2 and 176 
in N3) above a flux limit 0.135 mJy in the deeper observa- 
tions or 1.15 mJy over the shallower ones. We have selected 
sources for which we have optical data from the WFS (i.e., 
sources in Nl and N2). Multi-component sources (flagged 
in the catalogue as 'A', 'B' or 'C components) have been 
removed and only the calculated central position considered 
('T' in the catalogue). 
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Figure 4. Offsets between 15/im and optical positions (loft) and 
between radio and optical positions (right). 

4 OPTICAL IDENTIFICATIONS OF ELAIS 
SOURCES 

As shown in lLari et alJ ll200 j) . the positional errors in RA 
and DEC for the ELAIS sources result from the combina- 
tion of three quantities: the finite spatial sampling (cts), the 
reduction method {oy-) and the pointing accuracy (o"p). The 
simulation work carried out in the ELAIS fields yield next 
relations: 

a,+,(RA) = 1.00 + 17.17 X e(-°'5^^s/'~') (1) 
cr,+,(DEC) = 1.06 + 1.21 X fJ-'^^-^^xS/^) ^-2) 

These equations have been used to estimate the posi- 
tional errors due to mapping and reduction method as a 
function of each source's signal-to-noise ratio, S/N. 

The errors introduced by uncertainties in the ISOCAM 
pointing have been estimated by correlating t he ISO source s 
with the USNO catalogue of optical objects jMonetill993) . 
For each raster, both catalogues have been correlated us- 
ing a maximum search distance of 12". The median of the 
offsets values for all the ISO-USNO associations have been 
calculated,crp(RA), (Tp(DEC). Each source position has then 
been corrected for the offset found for each raster. The final 
positional error is then: 

oIa = 'Ts+r(RA) + a^(RA) -t- 0.4^ (3) 
'TDEC = ':r3V(DEC) + ai(DEC)-f0.4' (4) 

where a cr = 0.4" has been added to account for the optical 
errors. Typical positional errors are cr ~ 3". 

The correlation between ELAIS 15pim sources and op- 
tical obj ects has been carried out usin g a likelihood ratio 
method JSutherland fc Saunders^ll992^ ■ similar to the one 
which has been successfully applied to the identification of 
15/jm sources detected by ISO in the HDF-N by Mann ct al. 

The probability that an optical object of magnitude m 
is the true counterpart of a source with an error ellipse de- 
fined by its major axis, cti, and minor axis, 02, separated a 
distance r is given by 

^ Q{m)exp{-r^/2) 

27raicr2iV(m) ' 



where Q and N are the magnitude distributions of the 
sources and objects respectively. The reliability of such iden- 
tification is 

The identification process is carried out as follows. For each 
ELAIS 15^m source, all optical objects within a distance 
of 20" are selected. This list is our candidate list. For each 
object in our candidate list we calculate the values of the 
likelihood and reliability as given by equations above. The 
likelihood value (equation |5J gives us the probability that 
a candidate is the true optical counterpart of the source; 
but it only provides information about the probability of 
each candidate being the correct counterpart. The reliabil- 
ity (equation ISJ provides information about the number of 
candidates with high likelihood values. A candidate will have 
large values of likelihood and reliability if it is the only prob- 
able counterpart of a source. In case where there are multiple 
probable counterparts (in the sense of high likelihood) , they 
all will have low values of reliability. 

A candidate is selected to be the correct optical identi- 
fication of a ELAIS source when C > 0.8. Sources for which 
no candidates meet this requirement are fiagged as blank 
fields and represent ~ 8% of the total sample (section |7|l- 
Sources for which there are more than one candidate meet- 
ing this requirement, and have low values of reliability, are 
fiagged as having multiple possible counterpart and repre- 
sent ~ 8% of the sample. Finally bright stars, saturated in 
the WFS CCD data are also flagged; note that their astro- 
metric accuracy is poor. They represent ~ 22% of the sam- 
ple. Figure3](left) shows the offsets between ISO and optical 
identifications, excluding saturated stars. Uncertainties are 
well fitted by a Gaussian distribution of ct ~ 1". 

Optical identification of the ELAIS radio sources de- 
tected in the Nl and N2 areas is carried out using a similar 
procedure. In this case the radio catalogue provides measure- 
ments of the positional errors, so these are used in our like- 
lihood ratio algorithm. An optical counterpart is found for 
389 out of the 691 sources, i.e., there is a 44% of blank fields. 
Figure 01 (right) shows the offsets between radio sources and 
the optical counterpart. This provides a confirmation of the 
good accuracy of the astrometry of the 15 /im sources. 

Figure shows some example finding charts of sources 
detected at 15/im. These charts are 30" x30" in size and 
have been extracted from the r'-band images. Also shown 
a 3" radius circle centered on the ISO position. As shown 
in section |21 a large fraction of the sources are associated 
with bright galaxies. Example of sources with more than one 
plausible counterpart is given by e.g., J163636. 2-1-411049. 
The most likely counterpart is an object with r'=21.8 at a 
distance of 1.1 arcsec from the ISO position. Its likelihood 
is 0.988 and reliability is 0.773. There is another source of 
r'=23.0 separated 2 arcsec from the ISO position with like- 
lihood 0.960 and reliability 0.227. We choose the first one as 
the optical counterpart (based on higher values of likelihood 
and reliability) of the ISO source and flag it has having more 
than one plausible counterpart. Examples of these are also 
J160347.5+544552 and J160721.6-h535437. 

There are four ELAIS sources which have been merged 
in the final analysis 15 ^m catalogue (figure 6). Both possible 
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Figure 6. Double sources merged in the ELAIS 15/im catalogue. 
Mid-IR contours starting at 5(t and increasing at steps of Icr are 
shown plotted over the r'-band image. Each chart is 30" x30" in 
size. 



counterparts are given in the optical identification catalogue 
and the 15 /im flux is assigned to both. 

The table of optical identifications is available electron- 
ically in: http://www.ast.cam.ac.uk/~eglez/eid. Also in- 
cluded are postage stamps for all the sources in the r' band 
(grayscale and contours) as well as multiband grayscale find- 
ing charts. 

The format of the table is as follows: 

Column 1. International Astronomical Unit name of the 
source. Sources are listed in RA order. Sources detected in 
N2 are listed after those detected in Nl. 



Column 2. ISO coordinates (J2000) of the sources. 
Column 3. Coordinates (J2000) of the optical counterpart 
of the source. 

Columns 4 to 8. Aperture magnitude in U, g', r', i' and Z 
bands (aperture of radius 3.5 pixels - 1.16"). 

Column 9. Total r'-band magnitude as provided by SEx- 
tractor MAG_BEST parameter. 

Columns 10 to 15. Errors in previous magnitudes. 

Columns 16 to 21. Stellar classification as provided by 
WFS and SExtractor CLASS.STAR parameters. 

Columns 22, 23 and 24. Distance, Aa and /S.5 between the 
ELAIS source and the optical association. 

Columns 25 and 26. Likelihood of the identification for- 
mated as C/{1 -\- C) and reliability. 

Columns 27 and 28. Flux at 15/im and signal-to-noise ra- 
tio. 

Column 29. Optical fiag code as follows. Bl: source in a 
gap between chips or in the edge of a chip, B3: most likely 
optical identification when multiple counterparts, B4: blank 
field, B7: bright saturated star, B9: other plausible optical 
identification when multiple counterparts. 



5 MAGNITUDE DISTRIBUTIONS 

Figure |7] shows the magnitude distribution of the optical 
counterparts of ELAIS 15 sources. The first peak at 
r'~12-13 is caused by bright stars which are saturated in 
the WFS data. The second peak, due to extra-galactic ob- 
jects, is located at r'= 18. Most of the extragalactic sources 
are then associated with optically bright objects. A tail of 
faint objects is present at r'>20. 

In order to calculate the percentage of chance associa- 
tions, especially at the faintest magnitudes, we have simu- 
lated four catalogues in each Nl and N2 regions by offsetting 
20 arcsec in R.A. and DEC from the ISO position in four 
directions. The likelihood ratio technique was used to asso- 
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Figure 7. Magnitude distribution of 15 /xm sources in ELAIS Nl 
and N2. Hatched histogram shows the distribution of point-hke 
sources. Last bin, with right arrow symbols, represent the objects 
without optical counterpart. 
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Figure 8. Magnitude distribution of 1.4 GHz radio sources in 
ELAIS Nl and N2. Hatched histogram shows the distribution of 
point-like sources. Last bin, with right arrow symbols, represent 
the objects without optical counterpart. 



ciate the sources in these new catalogues in the same way as 
done for the real ones. The identifications obtained provides 
the distribution of chance associations. This percentage is 
~5% for r's;20 and increases to 20% at r'=24. 

Figure[7|also shows the magnitude distribution of point 
like objects. About 16% of the objects with magnitude r'>15 
are classified as point like. Their magnitude distribution has 
a peak at r'~19, a magnitude fainter than that for galaxies. 

The magnitude distribution of radio sources is shown 
in figure |S] Unlike the distribution of mid-IR sources the 
number of sources show an increase at fainter magnitudes. 
The number of point-like objects is very low (hatched his- 
togram) but the reliability of the CLASS_STAR parameter in 
SExtractor decreases at faint magnitudes and fails at r'^23. 



6 OPTICAL TO INFRARED FLUXES 

Figure 1^ shows the mid-IR to optical fluxes for the ELAIS 
sources. Stars have typically low mid-IR fluxes compared 
to their optical fluxes and are located in the region where 



their mid-IR flux is ten times smaller than their optical flux. 
Most of the extragalactic objects have mid-IR fluxes between 
1 and 100 times their optical fluxes. Acc ording to models 
of infr ared galaxies previously published iRowan- RobinsonI 
i200ll) . galaxies whose infrared emission is dominated by 
cirrus are located in the region < logSis/S'r' < 1. More 
infrared "active" galaxies, i.e., starbursts, AGN and Arp220- 
like objects are located in regions of mid-IR flux 10 times to 
100 times larger than their optical flux. 

Optical colours can be used to discriminate between 
AGNs and galaxies. As figure [THl shows AGNs are typically 
0.5 magnitudes bluer than galaxies (and as shown in 1101 
about 1 magnitude fainter). There is a population of point- 
like objects with galaxy-like colours, presumable highly ob- 
scure d. Their overall spectra can b e well fitted by a galaxy 
SED llRowan- Robinson et ahlboOSl) . 

The radio to optical flux ratio versus magnitude and 
versus optical colour are shown in figures 1111 and 1121 for 
the 1.4 GHz sources with optical counterpart. Those sources 
which show emission at 15/im are also marked with black 
symbols. Most of the 15 /im- radio coincidences are objects 
with magnitudes 15<r'<20 for which the cirrus component 
is the most plausible source of the infrared and radio emis- 
sion. 

Attending to their radio, mid-IR and optical properties 
(see figure I13II the bulk of the 15 ^m sources (those with 
< logSiz/ Sr' < 2) can be explained as a mixture of cirrus- 
dominated galaxies and starburst with smaller fractions of 
Arp220-like objects and AGN. 



6.1 Population with large infrared-to-optical flux 
ratio 

Figure ^ also shows the existence of a population of ob- 
jects with extreme mid-IR to optical fluxes and faint op- 
tical magnitudes. The nature of this population can orfly 
be studied with detailed spectroscopic foUowup observa- 
tions. Two of the point like objects in this region of the 
diagram have spectroscopic redshifts. The object associ- 
ated with ELAISC15_J164021. 5+413925 has a 15/.im flux 
of Sis — 2.98 mjy, an optical magnitude of r'= 20.3 and 
a redshift z — 0.60 (Perez- Fournon et al. 2003, in prep.). 
Its luminosity is L = 10^^ L©. The second object is asso- 
ciated with ELAISC15_J163655. 8-^405909. It has a 15/im 
flux of Sis = 0.72 mJy, an optical magnit ude of r'= 23.9 
and a redshift z = 2.61 iWiUott et alj|2003t) . Its luminosity 
is L = 10^2 Lq {Ho = 65, D.„, = 0.3, Ao = 0.7). 

Using infrared models of starbu rst galaxies and AGN , 
and adding reddening (as modeled bv lCalzetti et al.1 l)200n() 'l. 
we can explain the nature of this population. Objects with 
log 5*15 / Sr' > 2 can be explained as luminous starburst 
galaxies, with luminosities L ~ lO^'^ Lq, at redshifts z ~ 0.7 
and reddening A„ ~ 1 — 1.5. Fainter optical objects with 
log Sis/ Sr' ~ 3 have typically larger luminosities L > 
10^'^ 7/0, redshifts z ~ 1.3 and reddening Av ~ 2. Luminous 
AGN, with L ~ 10^^ Lq , also populate this area. 

The association of this sources with very reddened 
objects supports the assumption that the objects with 
Sis/Sr' > 100 in figure |n| may represent a new population 
of heavily obscured starbursts and type 2 AGN. 
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Figure 9. Optical to infrared flux ratio for 15^m sources with 
optical counterpart. Top figure shows the mid-IR flux versus op- 
tical magnitude. Point-like classified object are shown as black 
stars while galaxies are shown as gray circles. Lines of constant 
flux ratio log(Si5/5r/) = -1.0,0.0,-1-1.0,-1-2.0 are also shown. 
Blank flelds are displayed using a right arrow. Bottom flgure 
shows log(Si5/S'^/ ) versus optical magnitude. 



Table 1. ELAIS blank flolds with a radio source nearest than 3' 
from the ISO position. Fluxes are given in mjy. 



Name 


ISO Coords (J2000) 




'S'20cm 


J160734. 3+544216 
J163505. 4-1-412508 
J163511. 4+412255 


16 07 34.40 +54 42 15.6 
16 35 05.71 +41 25 11.2 
16 35 11.54 +41 22 57.4 


1.51+0.19 
0.83+0.09 
0.72+0.09 


0.37+0.02 
0.38+0.02 
0.94+0.02 



7 BLANK FIELDS 

A percentage of 8% of the 15/im sources (38 in Nl, 67 in 
N2) do not have an optical counterpart down to the opti- 
cal limits of data described in previous sections. All of them 
have infrared fluxes Sis < 3mJy, and only 3 show a plau- 
sible association with a radio source (table 0. Blank fields 
are probably the extreme version of the objects found at 
logSis/Sr' > 2. Using the same model galaxies as before, we 
find that starburst galaxies with luminosities L ~ 10^"^ L©, 
2 ~ 1. and Av ~ 3 — 4 may populate this region. The cause of 
the infrared emission and the starburst activity may be the 
merger of galaxies ir the formation of a protogalaxy although 
the accuracy of this hypotheses cannot be tested with these 
data. 
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Figure 10. Infrared to optical flux ratio versus optical colour 
g' — r' . Symbols as in figure 151 
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Figure 11. Radio to optical fiux ratio versus optical magnitude. 
Black symbols represent the sources detected also at 15/xm (dots: 
galaxies, stars: point-like objects). Hatched areas represent the 
range predicted by different infrared galaxy types in the range 
< 2 < 1.5 (redshift increases from bottom to top in each area; 
dashed line shows the location of 2 = 1.0). 
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Figure 12. Radio to optical flux ratio versus optical colour. Sym- 
bols as in figure im 
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Figure 13. Radio to mid-IR, flux ratio versus mid-IR to optical 
flux ratio. Symbols as in flgure 1111 Lines show the location of 
different infrared galaxy models.. 

8 SUMMARY 

The association of sources detected at 15/im in the ELAIS 
Nl and N2 areas with optical objects is presented. A 92% of 
the sample presents an optical identification to r'=24. The 
magnitude distribution presents a maximum at r'=18 and a 
tail which extends to fainter magnitudes. The distribution 
of point-like objects presents a maximum one magnitude 
fainter. The mid-IR to optical flux ratios, Si^/Sr', of the 
bright optical sources are in the range 1 to 10^ and can be 
explained using simple models of cirrus, starbursts, AGN 
and Arp220 spectral energy distributions. The tail of faint 
objects, show larger Si^/Sr', from 10^ to 10^ and can only 
be explained assuming large luminosities and obscurations. 
Point like objects show bluer g'-r' colour while higher mid- 
IR to optical flux ratio than galaxies. The remaining 8% of 
objects not identified are faint in the mid-IR , with a 15/im 
flux lower than 3 mjy. However, their mid-IR-to-optical flux 
ratio is larger than 10^ favouring the interpretation that they 
are associated with starbursts or AGNs at high redshifts and 
highly obscured. The identification of radio sources in the 
same areas is also presented. Their magnitude distribution 
shows an increase on the number of sources towards faint 
magnitudes. The number of unidentified objects is 44%. 
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